Geometrical terms in the effective Hamiltonian for rotor molecules 
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An analogy between asymmetric rotor molecules and anisotropic cosmology can be used to calcu- 
late new centrifugal distortion terms in the effective potential of asymmetric rotor molecules which 
have no internal 3-fold symmetry. The torsional potential picks up extra cos a and cos 2a contri- 
butions, which are comparable to corrections to the momentum terms in methanol and other rotor 
molecules with isotope replacements. 
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Geometrical ideas can often be used to find underlying order in complex systems. In this paper we shall examine 
some of the geometry associated with molecular systems with an internal rotational degree of freedom, and make use 
of a mathematical analogy between these rotor molecules and a class of anisotropic cosmological models to evaluate 
new torsional potential terms in the effective molecular Hamiltonian. 

An effective Hamiltonian can be constructed for any dynamical system in which the internal forces can be divided 
into a strong constraining forces and weaker non-constraining forces. The surface of constraint inherits a natural 
geometry induced by the kinetic energy functional. This geometry can be described in generalised coordinates by 
a metric gab ■ The energy levels of the corresponding quantum system divide into energy bands separated by energy 
gaps. A generalisation of the Born-oppenheimer approximation gives an effective Hamiltonian for an individual energy 
band. Working up to order , the effective quantum Hamiltonian for the reduced theory has a simple form [IHl] 

Hotf = ^igr'/vi^r w + V+ vgbo + vr, (1) 

where Pa = —ihd/dq'^, g"^ is the inverse of the metric, \g\ its determinant and V is the restriction of the potential of 
the original system to the constraint surface. Vgbo depends on forces orthogonal to the constraint surface. In the 
lowest energy band, the part depending on the normal mode frequency matrix u is given in Ref. [3], 

VaBo^^i^'^+Yf^t^i'^''^'^-'^'^^'^), (2) 

where Va is the Levy-Civita covariant derivative along the constraint surface. Terms depending on the anharmonicity 
of the potential can be found in Ref. [3 . The next term Vr is purely geometrical in nature [TJ [2] , 

Vn^^-h'R-^-h^K^ (3) 

depending on the intrinsic curvature scalar R and the extrinsic curvature scalar K of the constraint surface. The 
approximation to the effective Hamiltonian can be extended further whilst maintaining the symmetry under coordinate 
redefinitions by including additional momentum terms, such as 

n^tr{iy-^)R''''Wa'^b, n^ti{iy-^W''i^i^-^W''i^)VaVb... (4) 

The underlying geometry implies that the coefficients of these terms are universal. 

In molecular systems, constraining forces fix the length of the chemical bonds leaving the angular orientation of 
the molecule unconstrained. The potential V is the sum of the nuclear potential and the usual Born-Openheimer 
potential, V = Vnn + Vbo- Other terms in the effective Hamiltonian are caused by rotational- vibrationl coupling. The 
commonly adopted procedure is to write down a set of (non-geometrical) momentum terms and fit their coefficients, 
centrifugal distortion constants, using data from molecular spectroscopy [^. Many of these parameters can then be 
compared with their values calculated ab initio in terms of molecular structure constants [6] . The geometric expansion 
of the effective Hamiltonian ( 

Rotor molecules have additional internal rotational angles, as shown in figure 1. The spectra of rotor molecules 
are sensitive to the physical environment and to molecular structure, making these molecules important in many 
astrophysical and computational chemistry applications 7 . The aim of the present paper is to investigate the 
geometric potential term Vr in asymmetric rotor molecules. This term vanishes for molecules with an internal 3- 
fold symmetry, which are the only rotor molecules in which a formal theory of the centrifugal distortion has been 
thouroughly developed [5Hl2). Important aspects of this term are is its universality, and the fact that it can be 
calculated relatively simply from the molecular structure. The downside of the geometrical potential is that it is 
small compared to leading order terms in the potential which can only be evaluated ab initio by heavy-duty numerical 
calculations. 



2|4 ) provides a way to repackage this information in a coordinate independent way. 



The study of molecules with internal rotation goes back a long way, with many theoretical developments made 
in the 1950's (13) . Some recent methodology can be found in a review by Kleiner [7]. The basic set-up is shown 
schematically in figure [T] The molecule consists of two sets of atoms A and B which are free to rotate independently 
about a common molecular axis, with angle a between A and B. The molecule frame is centred on the centre of mass 
with z axis aligned parallel to the internal rotation axis |13j and the x axis is fixed with respect to B. The molecule 
coordinates are the Euler angles 9, (j) and ip of the molecular frame and the internal rotation a. The approach adopted 
here is to develop concepts which are, as far as possible, independent of the choice of coordinates. 



FIG. 1. Schematic diagram of a molecule with internal rotation. Atoms in A, the top, can rotate freely with respect to atoms 
in B, the frame. The dotted line denotes the internal rotation axis,. 



The mass the molecule is denoted by M and the mass of A by Ma- A vector a runs from the centre of mass to a 
chosen point C on the molecular axis so that the locations of the individual nuclei are given by r„ = R{d, 4>, + 
i?(0, 0, a„)a„), where rotation matrices R act on the vectors a„ fixed to A or i?, and a„ = a when n (z A, a„ = 
otherwise. The total moment of inertia of the molecule lij is computed in the molecule frame with origin at the centre 
of mass, whilst the moments of inertia of A and B are more conveniently centred on the rotation axis at C. 

The rotational kinetic energy of the molecule is [HI [15] 

T - ^I^.uj'io^ + p.uj'a + h^a^, (5) 

where a;' is the angular velocity in the body frame, and the components are given by 

- I^,, + /''.y - ^f(E, X cr) • (E, X cr) (6) 
p, = I^,, + M^(E, X a) ■ (E, X pa) (7) 
= I^.. - Ml (E, X pa) ■ (E, X pa)/M. (8) 

The vector pA runs from the point C on the rotation axis to the centre of mass of A. 

The kinetic energy defines an intrinsic metric g on the four dimensional configuration space. It is convenient to use 
the same notation for the differentials associated with the angular velocities, and w" = da for the internal rotation. 
The metric can be written in the form 

g = Jy-K + p*w")(cj^' + p'uj'') + N^Lu'^^. (9) 

Comparison with the kinetic energy shows that the vector = Pj and A^^ ^ ~ piPj, where is the matrix 
inverse of lij. The metric (|9| has a three dimensional symmetry group which can be classified under the Bianchi 
classification of Lie algebras as Bianchi type IX [T64 . This symmetry corresponds to rotating the molecular coordinate 
frame and allows us to change frame to another coordinate system whenever this is convenient. (This has nothing to 
do with the molecular symmetry group) . Bianchi IX metrics also arise in the study of cosmological models [TT^ . An 
anisotropic inertia tensor corresponds to an anisotropic cosmological model, and the internal angle is analogous to 
the cosmological time parameter, though a change from space to spacetime requires replacing iV^ by —N'^. 
The kinetic term in the effective Hamiltonian ([T]) is given by the Laplacian, 

- f = + 1/1-1/2^1/2 _ p^p^^ |J|l/2^1/2 _ ^p^^ (10) 

where F = l/2N^ and Pi are the usual angular momenta. The momenta can also be written as Pa = —ihDa, where 
Di are the derivatives defined by the duality relations oj'^{Dj) — and Da = 8/ da. The ordering of the derivatives 
in this way ensures that the quantum theory preserves the Bianchi IX symmetry. Note that whenever the top part 
part of the molecule has C„>2 symmetry about the rotation axis, then F and |/| and are all independent of the 



4 



internal angle and the factor ordering becomes unimportant. In this case it makes sense to choose a frame in which 
selected components of lij or vanish [13] . 

The geometrical potential terms in the effective Hamiltonian are invariant under the frame rotations up to relocation 
of the origin of the internal rotation angle. The curvature scalar may be calculated using general formulae for Bianchi- 
type metrics in Ryan and Shepley [ 17^ . If I denotes the inertia matrix, then 

R = -2F^/^D^lr {i-^F^/^dS) - \f {trl-^D^lf - ^Ftr (l-^D^l)^ (11) 

+2trl-i + ^trItrI-2 _ itrl(trl-i)^ (12) 

The remaining extrinsic curvature terms can be related to the tangent vectors = DaVn, a — i or a. The components 
of the intrinsic metric ^ are 

9ab = ^ mn{DaYn) ' (-Dhl"n) 
n 

The extrinsic curvature Kab is the normal projection of DbBa, 

Kab = DbDaTn " g""^ ^ mp{DbDaYp) ■ DjCp DdVn- (13) 

p 

Its square is given by 

7^2 ab cd ab cd ef /l /l\ 

K g Kabcd-g g g-'KabeUcdf- (14) 

where 

Kabcd = ^ m,i{DaDbr,i) ■ {DcDdVn), (15) 

n 

Kabc = ^ mn{DaDbVn) ' DcVn- (16) 

n 

It is possible to express both of these tensors entirely in terms of the inertia tensors 1^, pi and I^. If there is a C„>2 
symmetry about the rotation axis then the inertial tensors and the curvature terms are constant. The interesting 
cases are therefore ones in which the molecule has at most a C2 symmetry or in which the rotation axis is displaced 
from the axis of symmetry. 

The potential in the effective Hamiltonian has contributions from the zero-point vibrational energy and the extra 
vibrational term in ([2]), which uses the covariant derivatives of the normal modes 

Vii//j = {ui - vj) ^ m„i?(0, 0, an) d„/ x d„j. (17) 

n 

d„/ are the displacements of molecule n in normal mode / with frequency vi. This paper focusses on the less familiar 
geometrical terms, and the analysis given above shows that these have no dependence on the vibrational modes. 

For the most general type of asymmetric rotor, in an arbitrary molecular frame, the geometric potential has a 
Fourier series expansion 

C30 

Vr = (a„ cos na + 6n sin no) . 

The cofhcients can be combined to form a set which is indepenent of the frame, by taking Vg = ag and Vn = 
2(a^ + 6^)1/2. These reduce to the usual Fourrier series coefficients when 6„ = 0, which happens, for example, if a 
reflection symmetry has been used to align the origin of the molecular frame. 



Results for the geometric potentials Vr and the kinetic term F (see eq. (10)) of a selection of molecules with 
isotoptic replacements are presented in table 1. Breaking C3 symmetry results in geometric terms which are typically 
in the range 0.1-0.3 cm^^. Some of the molecules have non-vanishing potentials even though they posses C3 symmetry, 
due to misalignment of the rotation axis and the axis of symmetry. The components Vi and Fi split doublet into 
triplet spectral lines in the rotational spectra of molecules such as CD2HCHO [19]. The splitting due to the Fi 
(kinetic) term is known to be insufficient to explain the data without adjusting the alignment of the rotation axis |14| 
or introducing empirical centrifugal distortion terms 120) . 
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TABLE I. Fourrier series components of the geometric potential term Vr — jh^R— ^h?K^ and the kinetic function F for a 
selection of molecules with internal rotation. These use identical molecular structures for each isotopomer using data obtained 
from [18] . Q 

Isotopomer V\/hc{cm~^) V2/hc{cm~^) V-i/hc{cm~^) Fo/hc{crn^^) Fi/hcicm^^) F2/hc{cm~^) 



Ethylene 




0.000 


0.279 


0.000 


21.964 


0.000 


0.000 




CH2CD2 


0.000 


0.284 


0.000 


16.476 


0.000 


0.000 




CHDCHD 


0.126 


0.393 


0.003 


15.203 


0.496 


0.048 


Methanol 


CH3OH 


0.036 


0.001 


0.000 


27.639 


0.129 


0.000 




CH2DOH 


0.105 


0.154 


0.000 


26.392 


0.496 


0.060 




CH2D^*0H 


0.101 


0.137 


0.000 


26.216 


0.470 


0.051 




CHD2OH 


0.040 


0.110 


0.000 


25.480 


0.244 


0.050 




CHD2OD 


0.040 


0.126 


0.000 


15.187 


0.244 


0.049 


Peroxymethyl 


CH3O2 


0.000 


0.000 


0.000 


6.845 


0.000 


0.000 




CH2DO2 


0.043 


0.291 


0.001 


5.546 


0.359 


0.062 




CHD2O2 


0.015 


0.185 


0.001 


4.684 


0.232 


0.009 


Nitrosomethane 


CH3NO 


0.017 


0.039 


0.000 


7.680 


0.037 


0.007 




CH2DNO 


0.052 


0.297 


0.000 


6.390 


0.412 


0.082 




CHD2NO 


0.037 


0.238 


0.001 


5.562 


0.302 


0.047 


Acet aldehyde 


CH3CHO 


0.011 


0.001 


0.000 


7.744 


0.034 


0.000 




CH2DCHO 


0.064 


0.210 


0.001 


6.431 


0.408 


0.059 




CHD2CHO 


0.041 


0.140 


0.001 


5.589 


0.301 


0.023 



^ A program to calculate these terms can be downloaded from http://research.ncl.ac.uk/cosmology/publications.html 



TABLE IL Fourrier series components of the geometric potential term Vr — jh^R — ^h^K^ and the kinetic function F for a 
selection of halogenated organic molecules with internal rotation. 

Halogen substituent Vi/hc{cm~^) V2/hc{cmr^) V3/hc{cm~^) Fo/hc{cm~^) Fi/hc{cm~^) F2/hc{cm^^) 



Ethylene 


CF2CH2 


0.000 


0.275 


0.000 


11.555 


0.000 


0.000 




CHFCHF 


0.594 


1.941 


0.054 


5.291 


2.530 


1.055 


Methanol 


CHF2HO 


0.284 


0.601 


0.005 


22.353 


0.233 


0.440 




CH2FHO 


0.526 


1.090 


0.002 


23.790 


1.318 


0.440 




CH2FFO 


0.330 


1.573 


0.059 


4.106 


1.155 


0.746 


Peroxymethyl 


CH2FO2 


0.178 


1.378 


0.057 


2.776 


0.899 


0.587 




CHF2O2 


0.145 


0.472 


0.022 


1.665 


0.034 


0.181 


Nitrosomethane 


CH2FNO 


0.274 


1.542 


0.063 


3.572 


1.084 


0.736 




CHF2NO 


0.189 


0.689 


0.038 


2.425 


0.105 


0.398 




CF3NO 


0.019 


0.034 


0.000 


2.043 


0.050 


0.024 


Acet aldehyde 


CH2FCHO 


0.266 


1.128 


0.032 


3.637 


1.120 


0.559 




CHF2CHO 


0.135 


0.473 


0.017 


2.485 


0.057 


0.279 



With large asymmetry, for example with the halogenated organic molecules in table 2, the geometrical potenials 
can be as high as 1-2 cm~^. 

The usefuUness, or otherwise, of these geometrical terms is largely dependent on the accuracy of the calculations for 
much larger Born-Openhiemer and vibrational potential terms. The ab initio calculations for acetaldehyde in Refs. 
[3Tl[31] are said to have an accuracy of around 2 cm^^ for the coefficient V3. These calculations include the zero-point 
vibrational energy hti i//2. Note that, in addition to the zero point energy and the geometrical potential there should 
also be the extra vibrational term in the effective Haniiltonian [3] , 



(18) 
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Finally, there are numerical schemes available which go beyond the Born-Openheimer approximation by quantising 
the nuclei as well as the electrons (e.g. [23]). These methods give the energy of the ground state of the molecule, 
and by comparing different equilibium configurations for the internal rotation they can also give an effective barrier 
height, but they do not so far give information about the Fourier components of the effective potential. 
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